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Astrophysical analysis of the measurement of (a, 7) and (a,n) cross sections of 69 Tm 
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Reaction cross sections of 169 Tm(a, 7 ) 173 Lu and 169 Tm(a,n) 172 Lu have been measured in the 
energy range 12.6 < E a < 17.5 MeV and 11.5 < E a < 17.5 MeV, respectively, using the recently 
introduced method of combining activation with X-ray counting. Improved shielding allowed to 
measure the (a,y) to lower energy than previously possible. The combination of (0,7) and (a,n) 
data made it possible to study the energy dependence of the a width. While absolute value and 
energy dependence are perfectly reproduced by theory at the energies above 14 MeV, the observed 
change in energy dependence at energies below 14 MeV requires a modification of the predicted 
a width. Using an effective, energy-dependent, local optical a+nucleus potential it is possible to 
reproduce the data but the astrophysical rate is still not well constrained at 7-process temperatures. 
The additional uncertainty stemming from a possible modification of the compound formation cross 
section is discussed. Including the remaining uncertainties, the recommended range of astrophysical 
reaction rate values at 2 GK is higher than the previously used values by factors of 2 — 37. 

PACS numbers: 26.30.Ef 26.50,+x 24.60.Dr 25.40.Lw 25.55.-e 29.30.Kv 
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I. ASTROPHYSICAL MOTIVATION 

Two neutron capture processes, the s and r process, 
are required to produce the bulk of natural nuclides above 
Fe [THl]- These two processes cannot, however, create 35 
neutron-deficient, stable, rare isotopes between Se and 
Hg, which are termed p nuclei. Photodisintegration of 
stable nuclei in the O/Ne shell of massive stars during 
a core-collapse supernova explosion has been suggested 
as a production mechanism for these nuclei |^-|6|- Such 
a so-called 7 process commences by sequences of (7,n) 
reactions which are replaced by (7,p) and (7,0:) reactions 
when reaching sufficiently neutron-deficient nuclides in 
an isotopic chain Q- 

Two mass regions have remained problematic when ex- 
plaining the production of p nuclei by the 7 process in 
core-collapse supernovae: the lightest p nuclei with mass 
numbers A < 100 and those in an intermediate region at 
150 < A < 165 are underproduced @, H H-Q3. While 
the explanation of the light p nuclei most likely requires 
a different astrophysical model, the problem in the in- 
termediate mass region may still be solved by improved 
reaction rates. 

For the 7 process, photodisintegrations happen in the 
plasma temperature range 2.0 < T < 3.0 GK. The tem- 
perature is tightly constrained by the necessity to pho- 
todisintegrate the lighter, more tightly bound seed nu- 
clei while also retaining heavy p nuclides. Several lay- 
ers with slightly different temperatures contribute to p 
nucleosynthesis in a star. The rates for heavier nuclei 
have to be known in the lower part of the tempera- 
ture region because they would be destroyed completely 
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at slightly higher temperatures. Moreover, (7,0) reac- 
tions have been found important for the intermediate 
and heavy mass region, whereas (7, p) dominates in the 
lighter mass region of the 7-process [111 [T2I ] . The photo- 
disintegration reaction rates are usually computed from 
capture rates by applying the reciprocity principle of stel- 
lar rates [13l Il4j . Therefore a measurement of a capture 
may determine also the photodisintegration rates, pro- 
vided the g.s. contribution to the stellar rate is large. 
Due to the Coulomb barrier, the (a/y) reaction cross 
sections are tiny at astrophysical energy and thus cur- 
rently unmeasurable. Going to as low energy as possible, 
however, already some discrepancies between data and 
predictions have been found in previous measurements. 
Low-energy a capture on heavy nuclei seems to be of- 
ten overpredicted, although there is not yet enough data 
to draw a conclusive picture. It underlines, neverthe- 
less, that the underproduction of p-nuclides in the range 
150 < A < 165 observed in stellar models may have a 
nuclear physics cause. 

Investigations of rates for the 7 process are not only im- 
portant for nucleosynthesis in core-collapse supernovae. 
Simulations of the thermonuclear explosion of a White 
Dwarf (type la supernova) also found p nuclei being pro- 
duced in a 7 process [Hj]. Regardless of the site, 7 pro- 
cess studies require a sound determination of the relevant 
astrophysical reaction rates by nuclear physics investiga- 
tions. The 7-process reaction networks include hundreds 
of nuclei and thousands of reactions, mostly on unsta- 
ble nuclei. All of the reaction rates are predicted in 
the Hauser-Feshbach statistical model of nuclear reac- 
tions [H, Since (7,p) and (7,0:) reactions occur at 
unstable isotopes, current experimental techniques have 
to aim for testing reaction model predictions at stability 
and to provide the data for a global improvement of these 
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models and their input. 

The nucleus 169 Tm is not a p nuclide but it is close to 
the problematic mass range. Very few a-induced reac- 
tion data are known in this mass range and none close 
to the astrophysically relevant energy range. This made 
169 Tm an interesting target for investigation using the 
newly introduced method of activation with subsequent 
X-ray counting, supplementing the convential 7-counting 
methods. Details of the experimental method and first 
results were already published in [IH, [l||. Here we in- 
troduce the additional data, extending the (a,j) cross 
sections to lower energies, and focus on a discussion of 
the implications for constraining the astrophysical reac- 
tion rate for (a, 7) at 7-process temperatures. 



II. EXPERIMENTAL PROCEDURE AND 
RESULTS 

The first results on the studied 169 Tm+a reactions 
have been published already in abbreviated form [l8[ 
and a full description of the experimental technique has 
been presented in Recently, the shielding around 

the LEPS detector has been refined and thickened (de- 
tails can be found in [20l j ) this way the average labo- 
ratory background count rate is reduced to about 1.7 
x 1(T 3 1/ (keV s) at the 20 - 80 keV energy region. 
Using this improved shielding it was possible to study 
the 169 Tm(a,7) 173 Lu reaction at even lower energies, at 
E a = 13.0 and 12.6 MeV. The experimental approach 
was similar to the one published in [3, [l!| , here only ad- 
ditional information relevant for the recently measured 
cross sections is given. 

The thicknesses of the thulium targets - produced via 
vacuum evaporation onto 2 /im thick high purity Al back- 
ings - were 331 and 377 /J,g/cm 2 , corresponding to 1.18 
x 10 18 and 1.34 x 10 18 atom/cm 2 , respectively. The 
number of target atoms has been derived using weigh- 
ing and the PIXE method. A single irradiation using a 
13.0 MeV He ++ beam provided by the cyclotron accel- 
erator of ATOMKI was carried out, the total number of 
the He ++ particles impinging on the targets was 8.76 x 
10 17 . Similarly to [Il|,[22[ an Al energy degrader foil was 
placed between the two thulium targets, the degrader foil 
thickness was determined by measuring the energy loss 
of alpha particles emitted by a 241 Am radioactive source 
and by weighing. The energy losses in the thulium layer, 
in the backing and in the energy degrader foil were cal- 
culated using the SRIM code |23l |. 

After the end of the irradiation the activity of the two 
samples has been measured with the LEPS detector in far 
geometry to determine the 169 Tm(a,n) 172 Lu cross section 
by counting the yield of the emitted characteristic K Q i_2 
X-rays. The length of these countings was about a day 
and were repeated once more about 3-4 days later. The 
X-ray counting was carried out again in close geometry 
about 14-17 weeks later to determine the number of the 
173 Lu isotopes produced in the 169 Tm(a, 7) 173 Lu reac- 



TABLE I. Extended experimental S factors (and recently 
measured cross sections) of the 169 Tm(a,n) 172 Lu reaction; re- 
cently measured entries are indicated by bold numbers. 



E a -Ec.m. 

(MeV) (MeV) 



lw Tm(a,n) 1,z Lu 
(MeV) (10 28 MeV barn) 



"Tnifa,!!) 1 
recent data 
(/ibam) 



TuF 11.21 ±0.057 
11.85^ 11.55 ±0.059 
12. T 11.90 ±0.061 



12. 5 a 



12.19 ± 0.062 



11.153 
11.267 



11.491 
11.609 



11.839 
11.961 



12.128 
12.252 



12.6" 12.28 ±0.087 



12.193 
12.367 



13.0° 12.68 ±0.066 



12.614 
12.746 



13.5° 



14. a 



15. a 



15. 5 a 



16. a 



13.16 ± 0.066 



13.66 ± 0.069 



14.63 ± 0.075 



15.12 ± 0.077 



15.61 ± 0.079 



16. 5 a 16.10 ± 0.081 



17. a 16.59 ± 0.084 



17. 5 a 17.08 ± 0.086 



13.094 
13.226 



13.591 
13.729 



14.555 
14.705 



15.043 
15.197 



15.531 
15.689 



16.019 
16.181 



16.506 
16.674 



16.994 
17.166 



77.210 ± 8.127 
51.778 ± 5.450 



62.822 ± 6.746 
42.309 ±4.543 



54.544 ± 6.150 
36.906 ±4.161 



42.697 ± 3.799 
29.116 ± 2.591 



42.956 ± 4.024 
25.252 ± 2.366 



33.907 ± 3.218 
23.096 ± 2.192 



30.076 ± 1.984 
20.921 ± 1.380 



21.484 ± 1.638 
15.010 ± 1.144 



13.789 ± 1.087 
9.704 ± 0.765 



10.894 ± 0.892 
7.730 ± 0.633 



6.922 ± 0.546 
4.949 ± 0.390 



5.149 ± 0.376 
3.708 ± 0.271 



3.731 ± 0.276 
2.695 ± 0.199 



2.537 ± 0.188 
1.845 ± 0.137 



9.82 ± 0.92 



27.5 ± 2.61 



a taken from QJ, 
b this work 

c remeasured in this work; average between the remeasured value 
and the one from [lSl , weighted by the statistical uncertainties; 
supersedes [lcf 

d Average values weighted by the statistical uncertainties are 
given when two irradiations were carried out at the same 
energy. 



tion. During the cooling period of more than 14 weeks 
the 172 Lu activity of the targets decreased by a factor of 
about 16000, therefore the observed X-ray yield belongs 
solely to the decay of the 173 Lu. The length of these 
countings was about 3 weeks. 

The experimental data are shown in Tables U and [TTJ 
Consistent results were obtained for the cases of two ir- 
radiations at the same energy. At these energies the 
average cross section values weighted by the statistical 
uncertainty are given. Furthermore, to test the results 
based on X-ray counting, the activity of two samples 
were measured at the LNGS deep underground labora- 
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TABLE II. Extended experimental S factors (and recently 
measured cross sections) of the 169 Tm(a,7) 173 Lu reaction; re- 
cently measured entries are indicated by bold numbers. 



with 7/ being the Sommerfeld parameter 
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094 


/i n i ocn _L /in on 
401.850 ± 40.olz 








13 


226 






14.0^ 


13.66 ± 0.069 














13 


.591 


213.200 ± 18.109 








13 


.729 


148.950 ± 12.652 




15.0° 


14.63 ± 0.075 














14 


.555 


61.002 ± 5.995 








14 


705 


42.929 ± 4.219 




15. 5^ 


15.12 ± 0.077 














15 


043 


35.574 ± 3.655 








15 


197 


25.241 ± 2.593 




16.0 ad 


15.61 ± 0.079 














15 


.531 


18.394 ± 1.421 








15 


689 


13.152 ± 1.016 




16.5^ 


16.10 ± 0.081 














16. 


.019 


8.986 ± 0.694 








16. 


181 


6.471 ± 0.499 




17.(£ 


16.59 ± 0.084 














16. 


.506 


4.689 ± 0.376 








16. 


074 


3.387 ± 0.272 




17.5^ 


17.08 ± 0.086 














16. 


994 


2.679 ± 0.198 








17. 


166 


1.948 ± 0.144 





a taken from [T|, El 
b this work 

c remeasured in this work; average between the remeasured value 
and the one from [lgt | . weighted by the statistical uncertainties; 
supersedes [iH| 

d Average values weighted by the statistical uncertainties are 
given when two irradiations were carried out at the same 
energy. 

c counted at LNGS; the average (weighted by the statistical 
uncertainty and the uncertainty of the detector efficiency) of 
the measured values is given 



tory [19J. For these irradiations, the average weighted 
by the uncertainty of the measured cross section values 
is given here. The quoted uncertainty in the E c . m . val- 
ues corresponds to the energy stability of the beam and 
to the uncertainty of the energy loss in the target. The 
uncertainty of the cross section is the quadratic sum of 
the following partial errors: efficiency of the HPGe and 
LEPS detectors (6 and 4 %, respectively), number of tar- 
get atoms (4%), current measurement (3%), uncertainty 
of decay parameters (< 5 %) and counting statistics (0.5 
- 13%). The uncertainties given for the averaged values 
are the variances of the weighted means. 

In astrophysical investigations it is common to quote 
the astrophysical S factor. The cross section a(E) and 
the astrophysical S factor S(E) at cm. energy E c . m . are 
related by 



V 



2E r _, 



(2) 



S(E c . m .) — E c , m a(E c , m ,)e 



27tt; 



The charge numbers Z p , of projectile and target, re- 
spectively, and their reduced mass /j, enter the Sommer- 
feld parameter. Since the energy enters the calculation 
of the S factor also via the Sommerfeld parameter, the 
inclusion of the errors on cross sections and energies is 
not straightforward when computing the S factor from 
measured cross sections. Because of this, a range of S 
factors has to be given for each cm. energy Es in Tables 
U and HH evaluated at the lower and upper limit of the 
energy range defined by the errors on E c ra . The error 
bars on a energy and cross section translate into an error 
region for the S factor which is of trapezoid shape, with 
its four corners given by the upper and lower limit of the 
S factor at each energy E$- This error trapezoid is also 
shown in Figs. [4] and [5] 



III. DISCUSSION OF IMPLICATIONS FOR THE 
ASTROPHYSICAL REACTION RATE 

A. Relevant energy range and sensitivities 

In principle, reaction data can be used to test predic- 
tions of the cross sections and the resulting reaction rates 
for astrophysics in two ways. A direct comparison to data 
is useful if they reach the astrophysical energy range or, 
at least, are taken at energies where the cross sections 
show a similar sensitivity to the reaction widths as at 
astrophysically relevant energies. If such data are not 
available, the combined data of different reaction chan- 
nels may be used to extract information on the quantities 
determining the astrophysical reaction rates. 

The astrophysical energy window for 169 Tm(a,7) is 
6.3 - 9.2 MeV at 2 GK and 7.8 - 10.7 MeV at 3 GK 
[24l ]. The present data come close to the energy window 
but do not reach it. The reaction rate is only sensitive to 
the a width because it is the smallest width in the rele- 
vant energy region due to the Coulomb barrier (2f| . The 
sensitivity of the reaction cross section, shown in Fig. 
[TJ depends complicatedly on energy For a discussion of 
the sensitivity definition, see (T3, |25|. The sensitivity s 
is defined in such a way that |s| = 1 implies a change 
in the cross section by the same factor as the width was 
changed. When s < 1 the cross section values change 
opposite to the variation, i.e., they decrease when the 
corresponding width is increased. 

Below the neutron threshold, the (0,7) reaction is only 
sensitive to a change in the a width. Therefore it is sensi- 
tive to all uncertainties involving the optical a+nucleus 
potential. Above the threshold, it is still sensitive to 
the a width but also uncertainties in the neutron- and 
7-widths become increasingly important with increasing 
(1) energy. At E CAn , > 12.8 MeV, the sensitivity to these 
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FIG. 1. Sensitivity of the reaction 169 Tm(a,7) 173 Lu to a vari- 
ation of 7-, neutron-, and a-widths (variation of proton widths 
is not shown as the cross sections are insensitive to it) as 
function of cm. energy in MeV; the astrophysically relevant 
energy range 7.2 - 10.2 MeV at T = 2.5 GK is shown by the 
shaded region. 



FIG. 2. Sensitivity of the reaction 169 Tm(a,n) 172 Lu to a vari- 
ation of 7-, neutron-, and a- widths (variation of proton widths 
is not shown as the cross sections are insensitive to it) as 
function of cm. energy in MeV; the shown energy range was 
deliberately chosen to be the same as in Fig. [T]to facilitate a 
comparison. 



quantities exceeds the one for the a width. Sensitivities 
larger than one appear because a particles may be re- 
emitted in the 7 cascade when the compound nucleus 
de-excites 25]. Such emissions are only included ap- 
proximately in the version of the reaction code used and 
therefore larger uncertainties in the predictions may be 
expected at those energies. 

The energy range 12.19 < E c . m . < 17.17 MeV of the 
measured (01,7) cross sections covers a region where three 
different widths are important. The situation is fur- 
ther complicated by the fact that in this reaction with 
highly negative Q values in all channels, additional reac- 
tion channels appear, such as the cascade emission of a 
particles. In consequence, it is hard to disentangle the 
different contributions to test the astrophysically impor- 
tant quantity, the predicted a width. Cross comparison 
with the 169 Tm(a,n) 172 Lu data helps in this task. Figure 
[5] shows that the (a,n) reaction is mainly sensitive to the 
a width at the higher energies. Only close to the thresh- 
old, a stronger sensitivity to changes in the neutron- and 
7- widths appears. 

B. Comparison of theory to data 

Armed with the insights from the previous section, we 
can interpret the comparison between experiment and 
theory shown in Figs. [3] - [5] for the (a,n) and (0,7) 
reactions. The experimental S factors in these figures 
are compared to theoretical values obtained with vari- 
ous settings of the statistical model code SMARAGD, 
version 0.8.4s [TU, Hf|- It is important to mention 
that SMARAGD uses an improved routine to solve the 
Schrodinger equation, e.g., compared to the code used in 
[13, H3] ! leading to different results especially at energies 



close to the Coulomb barrier (see The calculated 

SMARAGD standard value (labeled 'std') uses the op- 
tical potential by [Hf. It perfectly reproduces the (a,n) 
data above 14 MeV, as seen in Fig. [3J For completeness, 
the S factors obtained with more recent global a+nucleus 
optical potentials from p9l - l33l | are also shown. They ex- 
hibit a very different energy dependence and cannot even 
describe the measured data above 14 MeV, with the ex- 
ception of the potential by ;33| which uses an energy- 
dependent parameterization approaching the parameters 
of [28j at high energies (see Sec. IIII C|) . Therefore, the 
following discussion of (a, 7) predictions and a possible 
extension to astrophysical energies focusses on the opti- 
cal potential of [28J and variations of the potential from 
[13]. 

Regarding the (a, 7) reaction, the standard theory val- 
ues exceed the experimental S factors by factors of 2.5 
and four at the upper and lower end, respectively, of the 
measured energy range (Fig. 3]). Apparently, not only is 
the absolute value not reproduced but also the energy 
dependence is different. Comparing the prediction to the 
(a,n) data in Figs. [3] and [5j however, shows that a mis- 
predicted a width cannot be the culprit for the deviation 
at energies above 14 MeV. As mentioned above, theory 
agrees with the experimental S factor values at these en- 
ergies. Below 14 MeV, on the other hand, a different 
energy dependence is seen in the (a,n) data when com- 
pared to theory. At first glance, this may be attributed to 
a problem in the prediction of the neutron- or 7- width be- 
cause the S factor is increasingly sensitive to these widths 
towards lower energy, as depicted in Fig.[5J In fact, either 
an increase in the 7 width or a decrease in the neutron 
width by a factor of 5 allows to reproduce the (a,n) data 
across the full energy range. 

It has to be realized, however, that the sensitivities 
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FIG. 3. (color online) Experimental astrophysical S factors 
(exp) of 169 Tm(a,n) 172 Lu as function of cm. energy are com- 
pared to values calculated with the global o+nucleus poten- 
tials by [H (std), [U (DGG), [H 12] (FR), [H (AA), and 
M (S) • 



to neutron and 7 widths are opposite in the (a,n) and 
(a, 7) reactions. In consequence, any attempt to remove 
the discrepancy of theory with the (a,n) data by modi- 
fying the neutron and/or 7 channel will result in shifting 
the calculated (01,7) S factors to higher values, further 
away from the experimental ones. Using the above mod- 
ification factor of 5, the new theory prediction for (0,7) 
would be a factor of 20 above the experimental values 
although the energy dependence is then reproduced well. 
Therefore, we have to proceed differently. The (a,n) data 
prove that the a width is described well above 14 MeV. 
The deviations below 14 MeV, nevertheless, may still be 
partly due to an incorrect energy dependence of the a 
width. From the (01,7) data, a factor of 2.5 between the- 
ory and experiment is extracted from the values above 14 
MeV. This factor must come from incorrectly predicted 
7- or neutron- widths. It is inconclusive with the given 
data which of the widths is incorrect and therefore we 
can compensate this factor by either increasing the neu- 
tron width or decreasing the 7 width. This will be of 
no consequence for the astrophysical rate which only de- 
pends on the a width. After having modified the neutron 
or 7 width for a given a width so that it reproduces the 
(a, 7) S'-factors above 14 MeV, remaining deviations from 
both (a, 7) and (a,n) data must be due to the a width. 
The investigations discussed in the following and regard- 
ing the a+nucleus optical potential have been based on 
a 7 width renormalized by a factor 0.5. All theoretical S 
factor curves shown in Figs. |H O and [6]- except for the 
'std' values - include the modified 7 width. 



C. Extension to astrophysical energies 

An energy-dependent optical a+nucleus potential, 
which transforms into the successful potential of [28| at 



FIG. 4. (color online) Experimental astrophysical S factors 
(exp) of 169 Tm(a,7) 173 Lu as function of cm. energy are com- 
pared to values calculated with different inputs; see text for 
further details. 
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FIG. 5. (color online) Same as Fig.Hbut for 169 Tm(a,n) 172 Lu 



projectile energies well above the Coulomb barrier but 
has a shallower imaginary part at low energy, was sug- 
gested in 33]. The potential uses the same real part 
and the same geometry as the potential of [28| but the 
strength of the volume imaginary part is parameterized 
with respect to the Coulomb barrier Eq 



W = 



25 



I + e (0.9£ c -Bc.m.)/ac 



MeV, 



(3) 



with ac — 2 MeV. At high energy E c . m . in the a channel, 
the depth W will assume the standard value from 28] . 

The S factor curves obtained with this potential are 
labeled 'S2' in Figs. IU [5] The energy dependence of the 
potential is too strong when applied to this reaction, as 
can be seen in Fig. [5j although it fared well in describing 
the 141 Pr(a,n) 144 Pm data of [Hj]. The only free param- 
eter is ac, determining the strength of the energy de- 
pendence. The curves labeled 'S4', 'S5', and 'S6' use ac 
values of 4, 5, and 6 MeV, respectively. All can describe 
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FIG. 6. (color online) Ratios rg x /r* td of the astrophysical 
169 Tm(a,7) 173 Lu reaction rates rg x calculated with different 
potentials Sa; (x = 4, 5,6) to the rate r* td calculated with the 
standard potential by [23] as function of plasma temperature 
T; the temperature range relevant in the 7 process is marked 
by the shaded region. 



the present 169 Tm(a,n) data well and simultaneously also 
the capture data, the latter being shown in Fig. SJ 

It should be noted that using the above approach, 
the optical potential effectively includes everything which 
may alter the a width at low energy. The a width can be 
affected by the included low-lying levels and the optical 
potential QUll- The level scheme of 169 Tm at low exci- 
tation energy is well established and can be assumed to be 
sufficiently known. This leaves the possibility of an incor- 
rect energy dependence of the optical potential. Another 
possibility, however, would be an alteration of the com- 
pound formation cross section by direct processes. The 
formation cross section is determined by the «Q ap width 
(i.e. the width for resonant a capture on the ground state 
(g.s.) of the target nucleus) and affects all compound re- 
action channels. An optical potential is usually derived 
from fits to elastic scattering data. Its imaginary part 
describes the flux into non-elastic channels but without 
distinguishing the nature of these channels. If such direct 
processes were implicitly included in the standard optical 
potential, the formation cross section would be overesti- 
mated because it is implicitly assumed in the reaction cal- 
culation that all a flux missing from the elastic channel 
will contribute to the formation of a compound nucleus. 
Therefore, the modification of the optical a+nucleus po- 
tential suggested above has to be viewed as an effective 
correction, ignoring the cause for the correction. If some 
sort of direct process was acting at low energy, explicit 
inclusion of such a process in the reaction calculation 
would remove the need to alter the optical potential. 

Assuming that Eq. is the correct function for de- 
scribing the energy dependence, the limits on the astro- 
physical reaction rates for 4 < clq < 6 MeV can be calcu- 
lated at plasma temperatures relevant for the 7 process. 
The astrophysical reaction rates resulting from the use 



FIG. 7. Comparison of g.s. contribution X (X) (taken 
from [25[) and stellar enhancement factor /sef (SEF) of 
169 Tm(a,7) 173 Lu as function of plasma temperature. The 
7-process temperature range is marked by the shaded area. 



of the three potentials are compared to the rate obtained 
with the standard potential in Fig. [JJJ The new rates are 
at least a factor of 5 below our standard rate obtained 
with the potential of [28j. At T = 2 GK, the ratios be- 
tween newly calculated rates and our standard prediction 
range from 0.06 to 0.136 when increasing ac from 4 to 6 
MeV. 

Table ILTTI gives the stellar reactivities Na {&v)* for the 
reactions 169 Tm(a,7) 173 Lu and 173 Lu,(7,a) 169 Tm com- 
puted with ac = 4, 5, 6 MeV at different plasma temper- 
atures. Table IIVI presents the coefficients for the usual 
REACLIB parameterization fl7| . Due to the negative 
reaction Q value, the parameter fits had to be treated 
specially. They were obtained by converting the calcu- 
lated stellar (a, 7) reactivity to the one for stellar (7,0;), 
using the reciprocity relations between forward and re- 
verse stellar reactivities [HI, This was then fitted 
and the parameters were converted back to the one for 
the (a, 7) reaction as shown in [X7j | . This ensures better 
numerical accuracy than fitting a reaction with negative 
Q value. 

As discussed above, it has to be cautioned, however, 
that the presented reactivities still carry large uncertain- 
ties as the energy-dependence of the a width at energies 
close to the Coulomb barrier is still not well understood. 
Moreover, strictly speaking the experiment only provides 
an indication of how the aQ ap width is modified. Stellar 
reactivities include transitions from and to excited states 
of the target nucleus, too. Figure [75] also shows the ratio 
of an "ao only" rate to our comparison standard. The 
" cko only" rate was obtained by using the potential mod- 
ification of Eq. ([3]) only for a transitions from and to 
the g.s. of 169 Tm while using the standard potential for 
all others. Incidentally, this yields an almost unchanged 
rate at 7-process temperatures. (The results for differ- 
ent values of ac are indistinguishable in the plot.) The 
shown behavior can be understood by realizing that the 



7 



TABLE III. Stellar reactivity Na for 169 Tm(a,7) 173 Lu as function of plasma temperature T obtained with three different 

values for ac - 4, 5, and 6 MeV - labeled by S4, S5, and S6, respectively. 



T 


S4 


S5 


S6 


S6 (ao only) 


(GK) 


(cm 3 s _1 mole _1 ) 


(cm 3 s _1 mole _1 ) 


(cm 3 s _1 mole _1 ) 


(cm 3 s _1 mole~ 


0.50 


1.009 x 10~ bz 


2.600 x 10~ b2 


4.557 x 10" 52 


4.600 x 10~ bi 


0.60 


2.600 x 10" 47 


6.301 x 10~ 47 


1.086 x 10 -46 


1.059 x 10 -45 


0.70 


4.942 x 10~ 43 


1.147 x 10~ 42 


1.923 x 10~ 42 


1.817 x 10~ 41 


0.80 


1.550 x 10 -39 


3.454 x 10~ 39 


5.644 x 10 -39 


5.183 x 10 -38 


0.90 


1.308 x 10~ 36 


2.807 x 10" 36 


4.488 x 10~ 36 


4.014 x 10~ 35 


1.00 


4.089 x 10 -34 


8.502 x 10 -34 


1.332 x 10~ 33 


1.162 x 10 -32 


1.50 


1.594 x 10~ 25 


2.902 x 10~ 25 


4.214 x 10" 25 


3.270 x 10" 24 


2.00 


2.522 x 10~ 20 


4.129 x 10 -20 


5.693 x 10~ 20 


3.920 x 10~ 19 


2.50 


9.190 x 10~ 17 


1.392 x 10~ 16 


1.838 x 10" 16 


1.122 x 10" 15 


3.00 


3.733 x 10 -14 


5.371 x 10 -14 


6.857 x 10~ 14 


3.721 x 10 -13 


3.50 


3.533 x 10~ 12 


4.911 x 10~ 12 


6.111 x 10" 12 


2.995 x 10" 11 


4.00 


1.157 x 10~ 10 


1.565 x 10 -10 


1.909 x 10~ 10 


8.594 x 10~ 10 


4.50 


1.707 x 10~ 9 


2.250 x 10~ 9 


2.697 x 10" 9 


1.124 x 10~ 8 


5.00 


1.385 x 10~ 8 


1.777 x 10~ 8 


2.093 x 10~ 8 


8.060 x 10~ 8 


6.00 


2.721 x 10~ 7 


3.283 x 10 -7 


3.725 x 10~ 7 


1.197 x 10~ 6 


7.00 


1.998 x 10~ 6 


2.271 x 10 -6 


2.481 x 10~ 6 


6.434 x 10~ 6 


8.00 


8.180 x 10~ 6 


8.874 x 10 -6 


9.406 x 10" 6 


1.967 x 10~ 5 


9.00 


2.184 x 10~ 5 


2.302 x 10 -5 


2.393 x 10~ 5 


4.217 x 10~ 5 


10.00 


4.117 x 10~ 5 


4.272 x 10~ 5 


4.391 x 10~ 5 


6.901 x 10~ 5 



TABLE IV. REACLIB parameters obtained from fitting the reactivities shown in Table ITTT1 



S4 



sr. 



S(> 



(Tig) 



(oyy) 



(tv 



S6 (a only) 
7) (7, a) 



ao -2.380584 x 1CT -2.143877 x 1CT 

ai -2.285393 x 10 1 0.000000 

a 2 -8.573430 x 10 1 

a 3 2.974330 x 10 2 

a 4 -2.842263 X 10 1 

a 5 1.740137 

a 6 -7.211434 x 10 1 -7.061434 x 10 1 



-2.525444 x 10 z -2.288737 x 10^ -2.134617 x 10^ -1.897910 x 10 z 

-2.285393 x 10 1 0.000000 -2.285393 x 10 1 0.000000 
-1.132236 x 10 2 -6.521460 x 10 1 

3.424164 x 10 2 2.499602 x 10 2 

-3.081515 x 10 1 -2.456620 x 10 1 

1.864776 1.437111 

-9.447688 x 10 1 -9.297688 x 10 1 -5.338706 x 10 1 -5.188706 x 10 1 



-1.521031 x 10^ -1.284324 x 10^ 
-2.285393 x 10 1 0.000000 
0.000000 
1.174915 x 10 2 
-1.591620 x 10 1 
8.700436 x 10 _1 
4.644152 6.144152 



temperature dependence of the altered g.s. rate is folded 
with the contribution X of the g.s. transitions to the total 
stellar rate. Figure [7] shows values for the g.s. contribu- 
tion X as function of temperature, taken from [25j . It 
can be seen that X is small at 7-process temperatures. 
Often, the stellar enhancement factor (SEF) 

(av)* 

/SEF - 7 TgTs- (4) 

is quoted as a comparison between stellar reactivity and 
the reactivity computed from g.s. transitions alone. The 
SEF is unsuited, however, to judge the contribution of the 
laboratory cross section to the stellar rate, as discussed in 
[34| . The reaction 169 Tm(a,7) 173 Lu is a good example 
for that. As found in Fig. [Jj /sef ~ 1 in the relevant 
temperature range whereas X is much lower. 

The stellar reactivities and the REACLIB parameters 
for the "ao only" case and ac = 6 MeV are also given 
in Tables IIIII and IIV1 respectively. Without a further 
understanding of the nature of the a width modification, 
the upper limit of the stellar reactivity is given by the 
values of model S6 (ao only) in the last column of Table 



IIIII while the lower limit is provided by the reactivities 
for S4 (a c = 4 MeV). 



IV. SUMMARY AND CONCLUSION 

We have determined the reaction cross sections 
and astrophysical S factors of 169 Tm(a,n) 172 Lu and 
169 Tm(a,7) 173 Lu at low energies, using our newly devel- 
oped method of combining activation and X-ray count- 
ing. An improved shielding around the LEPS detector 
enabled us to measure the (a, 7) reaction down to 12.6 
MeV a energy, lower than before. 

The impact of the new results on the determination of 
the astrophysical reaction rates for a capture on 169 Tm 
were discussed, using the full range of data. The combi- 
nation of (a,n) and (a,j) data was essential to disentan- 
gle errors in the predicted a width from those in other 
widths appearing in the reaction. It was found that the 
energy dependence of the a width below 14 MeV is dif- 
ferent from what is expected by standard predictions. A 
modified, energy-dependent, local optical a+nucleus po- 
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tential was presented, able to describe both (a,n) and 
(a, 7) S factors well across the measured energy range. 

Using the local potential, stellar (0,7) and (7, a) reac- 
tivities were calculated. They were found to be consider- 
ably lower at astrophysical 7-process temperatures than 
predictions using a standard optical potential. More spe- 
cifically, they were factors of 7.4 — 16.7 below the rate 
calculated with the potential by [28]. Ambiguities in the 
extrapolation to low energies, however, require an uncer- 
tainty of a factor of 2 — 3 in the predicted rate, even when 
assuming that the shape of the energy dependence is un- 
derstood. Further (0,7) measurements below the neutron 
threshold would be able to reduce this uncertainty. 

A further problem remains with identifying the na- 
ture of the potential modification. This leads to an even 
larger uncertainty in the calculated rate. Again, more 
low-energy data and an extended database including a 
wider range of nuclei would help to shed light on this 
issue. In conclusion, the recommended reactivity at a 
7-process temperature of 2 GK is 0.06 — 0.95 times the 
SMARAGD reactivity using the potential by [Hj], which 
translates to 2.3 — 37.0 times the widely used standard 
values of [27| . thereby leading to an enhancement in the 
(7, a) rate with respect to the values given in (l7l |27j|. 

The case studied here is a good example for the restric- 



tions and possible pitfalls that can be encountered when 
deriving astrophysical reaction rates from experimental 
data. It not only shows the importance of further mea- 
surements of reaction cross sections involving low-energy 
a particles to allow global studies of suitable optical po- 
tentials but also the significance of properly accounting 
for both, changing width sensitivities as well as thermally 
excited states, when interpreting the impact on an astro- 
physical reaction rate. 
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